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Suzaku Observation of Abell 1555 and Abell 1558: Searching for 
Non-thermal Emission from Large Scale Structure Formation 



(N 

o 

(N 

Oh- 
< 
00 



o 

u 

43 

9 L,: 
6 

CO 

03 



> 
CO 
On 
CO 

O 

CN 



1 



Ryu Makiya, Tomonori Totani 
Dept. of Astronomy, KyotoUniversity, Kitashirakawa-Oiwake-cho, Sakyo-ku, Kyoto 606-8502 
makiya@kusastro.kyoto-u.ac.jp, totani@kusastro.kyoto-u.ac.jp 
Kazuhiro NAKAZAWA 
Department of Physics, The University of Tokyo, 7-3-lHongo, Bunkyo-ku, Tokyo 113-0033 
nakazawa@amalthea.phys. s. u-tokyo. ac.jp 

(Received 2000 December 31; accepted 2001 January 1) 
Abstract 

We report X-ray observations of two galaxy clusters Abell 1555 and Abell 1558 with Suzaku, which 
are included in a large scale filamentary structure and a supercluster, to search for non-thermal emission 
driven by shocks produced in structure formation. These two clusters are detected by Suzaku/XIS for the 
first time in the X-ray band of 0.5-7 keV. No significant flux is detected by HXD in the energy band of 
13-40 keV, and upper limits are reported. From the analysis of the XIS data, we find that the spectrum of 
A1555 is fit by a thermal plus power-law model, significantly better than a singlc-tcmpcraturc pure thermal 
spectrum. If this power-law component is due to inverse-Compton scattering, the fraction of total baryon 
energy imparted to non-thermal electrons is consistent with the typical value inferred from the observation 
of other clusters. However, other scenarios (e.g., under lying AGNs, multi-temperature thermal models) 
cannot be excluded and further investigation of this system is desired. Basic physical properties of A1555 
(e.g., total mass) are also reported. 
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1. Introduction 

Galaxy clusters are the largest gravitationally bound 
structures in the Universe. It contains ~ 10 15 M Q hot 
gas, galaxies, dark matter, and non-thermal components 
such as magnetic fields and cosmic-rays. Although the 
detailed properties of galaxy clusters are intrinsically in- 
teresting topics, they are also important to understand 
the formation and evolution of the large scale structure, 
and acceleration mechanisms of cosmic-rays. 

According to the cold dark matter scenario, which is the 
most successful theory of the structure formation, struc- 
tures grow up hierarchically from small scale to large scale. 
When an object collapses gravitationally and virializes, 
shock waves occurs and a fraction of the baryonic mat- 
ters in objects are accelerated into relativistic energies, 
if the shock waves are strong enough. Accelerated elec- 
trons scatter the Cosmic Microwave Background (CMB) 
photons up to hard X-ray and gamma-ray regions via 
the Inverse-Compton (IC) mechanism (for a recent review 
of gamma-rays from structure formation, see Blasi et al. 
2007). 

Although the IC emission due to the structure forma- 
tion has not been detected at high significance, the exis- 
tence of relativistic electrons in several clusters has been 
revealed by the observations of diffuse radio emission from 
galaxy clusters (Willson 1970; Giovannini & Feretti 2000; 
Govoni & Feretti 2004; Feretti 2005). These extended ra- 
dio emissions are called as "radio halos" , which is located 
at the cluster center, or "radio relics" at the cluster periph- 



ery. They have steep spectra and a cutoff at a few GHz, 
implying that the relativistic electrons account for the ra- 
dio emission. The radio halos or relics are found in the 
merging clusters, and hence these electrons are considered 
to be accelerated in the cosmological shock waves. 

Since cluster X-ray emission is dominated by a thermal 
component below ~ 20 keV and expected non-thermal 
emission has a low flux, observing the IC emission is a 
difficult task. Although great effort has been devoted to 
search for the non-thermal IC emission, results are still 
controversial. For example, Beppo SAX reported the de- 
tection of hard X-ray emission from Coma cluster (Fusco- 
Fcmiano et al. 1999; Fusco-Fcmiano ct al. 2004), but these 
detections are still a matter of debate (Rossetti & Molendi 
2004; Fusco-Femiano et al. 2007). Observations of galaxy 
clusters in gamma-ray band have also been performed by 
ground-based air Cherenkov telescopes and Fermi (e.g., 
Aharonian et al. 2009a; Aharonian ct al. 2009b; Aleksic 
et al. 2010; Ackermann ct al. 2010), but no significant de- 
tection of non-thermal emission has been achieved, only 
resulting in flux upper limits. 

Here we report X-ray observations of two galaxy clus- 
ters, Abell 1555 and 1558 by Suzaku, whose redshifts are 
z = 0.127 and 0.116 (Colafrancesco 2002). The original 
motivation of this observation was to examine a hypothe- 
sis that an unidentified EGRET gamma-ray source 3EGJ 
1234-1318, which is located at the region of large scale fil- 
amentary structure including these two clusters (Einasto 
et al. 2001; Kawasaki and Totani 2002; sec fig.l), is IC 
emission from electrons accelerated by structure forma- 
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tion shocks. Recently, 3EGJ 1234-1318 was detected by 
the new generation gamma-ray telescope Fermi as 1FGL 
J1231. 1 — 1410, but its location has changed outside of 
the region of our Suzaku observation, and furthermore, 
its origin has been determined to be a gamma-ray pulsar 
(Ransom et al. 2011; Maeda et al. 2011). 

However, the region around A1555 and A1558 is still 
interesting as a general study for a supercluster region 
searching for non-thermal emission clusters. There have 
been no targeted deep X-ray observations for these clus- 
ters. We will present the result of our search for non- 
thermal nature of X-ray emission, as well as a study of gen- 
eral physical properties of these two clusters. The Hard 
X-ray Detector (HXD; Kokubun et al. 2007; Takahashi et 
al. 2007) onboard Suzaku (Mitsuda et al. 2007) is charac- 
terized by its low detector background and wide field of 
view for hard X-rays at 10-50 keV (HXD-PIN) and 50- 
600 keV(HXD-GSO). Suzaku X-ray CCD cameras (XISs; 
Koyama et al. 2007) has a large effective area and a stable 
and low detector background in 0.5-10 keV. These detec- 
tors are suitable for observing diffuse emission from galaxy 
clusters. 

In section 2 we present observation logs and data re- 
duction method, and section 3 is devoted for describing 
the result of data analysis. In section 4 we discuss the 
non-thermal nature of this supercluster. We will sum- 
marize our results in section 5. Throughout this pa- 
per, the cosmological parameters of Qq = 0.3, Qa = 0.7, 
h = Hq/ (100 Mpc -1 km s _1 ) = 0.7, and the baryon density 
parameter Hb — 0.015/i -2 are assumed. Unless otherwise 
noted, all errors are at the 90% confidence level. 

2. Observation and Data reduction 

The observations were carried out on 2006 December 
12-13 and 2007 December 11-12 for four pointings. 
Observation pointings are located along the filamentary 
structure and covers two galaxy clusters, A1555 and 
A1558. Hereafter we call these four observation areas as 
region 1, 2, 3, and 4, respectively (see figure 1 and figure 
2). Details of the observation dates and exposure times 
are summarized in table 1. In latter two observations, the 
trouble in onboard data processing of XIS-0 occurred and 
no usable data was output for the half area, and therefore 
we do not use XIS-0 data in this study. 

The data for region 1 and 2 were processed with the 
Suzaku pipeline processing of version 2.0.6.13, and region 
3 and 4 were processed with the version 2.2.8.20. We used 
the calibration data files of 20081009, and HEASOFT ver- 
sion of 6.6.3. XISs data obtained with 3x3 and 5x5 edit 
modes for each observation and we combined them into 
one and applied standard filters as follows. Events with 
a GRADE of 0, 2, 3, 4, 6, and STATUS < 1024 were 
extracted. We selected good time intervals by removing 
the time that space craft located at the South Atlantic 
Anomaly (SAA) within 436 sec, the cut-off rigidity < 6.0 
GV, elevation angle from the earth rim < 5°, and the 
sun- lit earth rim < 20°. 

Major background in Suzaku observation is Non X- 



ray Background (NXB) and Cosmic X-ray Background 
(CXB). NXB images and spectra of the XISs were created 
using the ftool "xisnxbgen" (Tawa et al. 2008). Tawa et 
al. (2008) studied the NXB error on a typical observation 
lasting for a few days. At the 90% confidence level, it 
is 6.7% and 12.5% for the sum of two Front Illuminated 
(FI) XISs (XIS-0,3) data and the Back Illuminated (BI) 
XIS (XIS-1) data, respectively. In the following analysis, 
we utilize this value as a systematic error of NXB model. 
The CXB fluctuation can be modeled as ccxb/-^cxb C£ 
nj°- 5 S°- 25 (Condon 1974). Here, Q e is the effective solid 
angle and S c is the upper cut-off flux. From the HEAO-1 
A2 results, Shafer (1983) derived <t C xb/Icxb = 2.8%(1ct) 
with il e = 15.8 deg 2 and S c = 8x 10" 11 erg s" 1 cm" 2 . 
By scaling this result with XIS parameters, we obtain 
O e = 0.09 deg 2 and S c ~ 4.94 x 10- 14 erg s _1 cm" 2 , which 
is the lowest flux of point sources detected in this observa- 
tion. CXB fluctuation in XISs were estimated to be 9.65% 
(90% confidence limit). 

The HXD data were also processed in a standard way. 
We screened the data with following criteria: the cut-off 
rigidity is larger than 6.0 GV, the elapsed time after the 
passage of SAA is more than 500 sec and the time before 
entering the SAA is more than 180 sec, elevation angle 
from the earth rim > 5°. We used a public NXB model 
provided by the HXD team for the NXB of the HXD- 
PIN. The version of the model is "METHOD=LCFITDT" 
or "tuned" (Fukazawa et al. 2009). 

According to the Fukazawa et al. (2009), the NXB 
model reproducibility of blank sky observations separated 
into 10 ks exposures give distribution of 5.7% at the 
90% confidence level, including the contribution from the 
statistical error of typically 3.3% or larger and the ef- 
fect of CXB fluctuation (1.3% of the total background). 
Thus the systematic error of NXB model is calculated 
to be 4.5% at the 90% confidence level. We use the 
same model as Nakazawa et al. (2009) for the CXB 
of HXD; N(E) = 8.7 x lO" 4 ^ 1 ' 29 x exp(-£/40.0) in 
photons -1 s _1 keV -1 FoV -1 , where E is the photon en- 
ergy in keV. Details of the CXB model estimation are 
described in appendix 2 of Nakazawa et al. (2009). Using 
the same model of the CXB fluctuation with the XIS, 
we obtain the CXB fluctuation in HXD as 18%, with 
O e = 0.32 deg 2 and conservative upper cut-off flux of 
5 c ~8x 10" 12 erg s -1 cm -2 in 10-40 keV band. 

3. Results 

3.1. HXD results 

HXD-PIN detected no significant signals above the 
background in all regions and we only set upper limits 
in each observation region using 13 keV - 40 keV data. 
We converted residual signals (data-(NXB+CXB)) into a 
energy flux in 13 keV - 40 keV by assuming a power law 
emission with photon index T = 2.0 located at a center of 
HXD FoV. Results are summarized in table 2. Systematic 
errors are the sum of CXB and NXB model uncertainty as 
described in section 2. In all regions the residual signal is 
negative. It means that CXB and NXB models are slightly 
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Fig. 1. The cluster "richness" map for our observation region overlaid (blue) by X-ray photon count contour (white solid line). We 
use the blue passband data of the APM galaxy catalog to make the "richness" map. See Kawasaki and Totani (2002) for detail. 
X-ray photon count contour is obtained by XIS-3 in the energy range of 0.5 - 7.0 keV. The contour levels are 2, 2.4, 2.8, 3.4, 4.0 
(in 10~ 7 cts s _1 pix — 1 ). Solid and dotted green boxes are XIS and HXD FoVs, respectively. Plus signs are galaxy cluster detected 
by Kawasaki and Totani (2002). Dotted yellow circle represents the error circle of 3EG J1234-1318 (The location of this source has 
been changed by Fermi observation; see text). 




Fig. 2. XIS-3 0.5 keV - 8.0 keV mosaic image. Green boxes are XIS FoV. Images were corrected for exposure time and vignetting 
effect after subtracting NXB, and smoothed by a Gaussian kernel with a = 0.5 arcmin. Small green circles show the detected sources. 
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Table 1. Observation log. 


Obs. region 


Obs. ID 


Start (UT) 


End (UT) 


Exp. (ks)* 


Region 1 
Region 2 
Region 3 
Region 4 


801032010 
801031010 
802004010 
802005010 


Dec 12, 2006 22:16 
Dec 12, 2006 08:56 
Dec 11, 2007 01:27 
Dec 11, 2007 20:31 


Dec 13, 2006 18:38 
Dec 12, 2006 22:15 
Dec 11, 2007 20:30 
Dec 12, 2007 11:10 


25.0/18.7 
10.7/26.1 
23.2/26.3 
13.8/20.0 



* The XIS/HXD Exposure time after screening 









Flux [erg cm 2 s x ] 




Upper limit 


Rej 


jion 1 


(" 


-7.4±2.3±5.9) x 10" 


-12 


8.2 x 10" 12 


Rej 


;ion 2 


(" 


-0.7±2.8±6.0) x 10" 


-12 


8.8 x 10~ 12 


Rej 


;ion 3 


(" 


-0.1 ±2.4 ±6.0) x 10- 


-12 


8.4 x 10~ 12 


Rej 


;ion 4 


( 


-0.2±2.8±6.1) x 10" 


-12 


8.9 x 10~ 12 



Table 2. Measured flux and upper limits in 13 KeV - 40 
keV from the HXD data. Spectral shape is assumed to be a 
power law with photon index 2. The former error is statistical, 
and the latter is systematic, due to CXB and NXB model 
uncertainty. All errors are at the 90% confidence level. 

overestimated, and we derived the flux upper limits from 
statistical and systematic errors, assuming physical flux to 
be zero. The comparison of these limits with the expected 
non-thermal IC emission is discussed in section 4.1. 

We also check the HXD-GSO spectrum and confirmed 
that it is consistent with the NXB model. Thus we did not 
use them. 

3.2. XIS results 

3.2.1. Image analysis 

In figure 2, we show the XIS3 mosaic image. Images 
were corrected for exposure time and vignetting effect af- 
ter subtracting NXB, and smoothed by a Gaussian kernel 
with a = 0.5 arcmin. Diffuse emissions from A1555 and 
A1558 are clearly detected, which are the first detection 
in the X-ray band. We also selected seven circular re- 
gions centered on significant signals selected by eyes (rep- 
resented by green circles in figure 2). On the other hand, 
there arc no apparent emission associated with the fila- 
mentary structure (see figurel). 

The non-thermal diffuse emission is expected to be spa- 
tially extended. To examine whether the signals except 
for the A1555 and A1558 are diffuse or not, we compare 
the spatial expanse of these signals with the Suzaku point 
spread function (PSF) . Since the PSF changes along with 
the source location on the XISs FoV, we estimated the 
PSF at the each source position using the ftool "xissim" . 
We simulated the image of point-like sources with the 
same exposure time as our observation, and utilize the 
spatial photon distribution of simulated point source as 
the PSF. 

We fitted the PSF and photon distribution of observed 
source with a Gaussian function. If the dispersions of the 
PSF and the photon distribution of an observed source 
agreed within ler, we regard this source as a point like 
source, and otherwise we regard it as a diffuse source. 



As a result, one source is determined to be diffuse (named 
"SUZAKU J1231-1307" , see figure 2) and others are point- 
like sources. We present detailed spectral analysis of this 
diffuse source, as well as A1555 and A1558, in next sub- 
section. 

Physical properties of all sources such as positions, pho- 
ton flux, point-like or diffuse, and radio counterpart (if ex- 
ists) are summarized in table 3. To derive a photon flux 
of a source, we assumed a power law model. Since fitting 
calculation is unstable when photon index is free, pho- 
ton index is fixed at 2. We selected sources from NLAO 
VLA Sky Survey (NVSS) catalog (Condon ct al. 1998) 
which are located within 0.3 arcmin from the source cen- 
ter as a radio counter part, where 0.3 arcmin is the typ- 
ical position determination accuracy of Suzaku. We find 
that SUZAKU J1231-1325 and SUZAKU J1233-1328have 
a radio counterpart, NVSS J123136-132502 and NVSS 
J123258-132727, respectively. 
3.2.2. Spectral Analysis: I. Data Reduction 

We performed spectral fitting of three diffuse sources, 
A1555, A1558 and SUZAKU J1231-1307 to search for the 
non-thermal diffuse emission. For the spectral analysis 
of XIS data, rmf and arf files were created using Suzaku 
ftool "xisrmfgen" and "xissimarfgen" , respectively. Two 
XIS-FI (XIS0, XIS3) spectra and response files for each 
sources are summed up, and XIS-FI and XIS-BI (XIS1) 
spectra are fitted simultaneously. Energy band used in 
the spectral analysis is 0.5 - 7.0 keV. 

The background spectrum of each FoV were modelled 
as follows. First, we extracted the spectrum from the 
whole FoV and subtracting the sources and NXB, and 
then performed spectral fitting of this data with the CXB 
model. For the CXB model, we used a model including a 
power-law component with a fixed photon index T = 1.4, a 
fixed temperature (0.08 keV) thermal component, and an- 
other thermal component whose temperature is a fitting 
parameter around 0.3 keV (the so-called "Galactic compo- 
nents"). The power-law component suffers from Galactic 
absorption, and it was modeled by the wabs code. For the 
Ah i column density we use the value of Kalberla et al. 
(2005), which are 3.47, 3.41, 3.41 and 3.43 x 10 20 cm- 2 
for the regions 1, 2, 3 and 4, respectively. 

Results of spectral fittings to CXB are summarized in 
table 4, where the values of fitting model parameters (the 
three normalization parameters and a fitting temperature) 
are shown. In each region, our CXB model describes 
the observed data well (see the reduced x 2 in the ta- 
ble), and the 2-10 keV fluxes of the power-law compo- 
nent are consistent within 90% statistic and systematic 
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Table 3. The list of sources which are detected in this observation. See the text for detail. 



Name 


RA 


DEC 


Flux * 


Point/Diffuse 


Counter part 


SUZAKU J1232-1259 


12:32:18 


-12:59:41 


1.4 ±0.20 


Point 




SUZAKU J1232-1302 


12:32:00 


-13:02:55 


0.62 ±0.13 


Point 


- 


SUZAKU J1232-1304 


12:32:07 


-13:04:42 


0.49 ±0.12 


Point 




SUZAKU J1231-1307 


12:31:44 


-13:07:08 




Diffuse 




SUZAKU J1232-1317 


12:32:17 


-13:17:02 


1.0±0.19 


Point 




SUZAKU J1231-1323 (A1555) 


12:31:60 


-13:23:33 




Diffuse 




SUZAKU J1231-1325 


12:31:30 


-13:25:00 


0.78 ±0.28 


Point 


NVSS J123136-132502 


SUZAKU J1233-1328 


12:33:02 


-13:28:46 


0.65 ±0.24 


Point 


NVSS J123258-132727 


SUZAKU J1234-1334 (A1558) 


12:34:00 


-13:34:32 




Diffuse 





* Energy flux in 10 13 erg s 1 cm 2 , where errors are at the 90% confidence limit. Flux of three diffuse sources are shown in table 5. 



error with the best- fit value derived using ASCA, which 
is 5.67 ±0.04 x 10 -8 erg cm" 2 s" 1 str" 1 (Kushino ct al. 
2002) . It has been known that solar- wind charge exchange 
with H I in the Earth's geocorona accounts for a part of X- 
ray background below IkeV (e.g., Cravens 2000; Fujimoto 
et al. 2007). We do not consider detailed characterization 
of this component in this paper since our CXB model well 
describes the spectra of background regions. (Note that 
this component may be implicitly included in the model- 
ing of the Galactic component.) 

Then, we simulated the CXB spectrum to produce X- 
ray counts using the obtained parameters of the CXB and 
arf files at each source position. Finally, we combined the 
simulated CXB and NXB spectrum, and utilized it as a 
background spectrum in the spectral fitting of sources. 
3.2.3. Spectral Analysis II. Model Fittings 

We fitted the spectra of the three diffuse sources by 
thermal, power-law, and thermal plus power-law mod- 
els. For the thermal model we use the apec code in the 
XSPEC (vll.3.2). For the value of 7V H i column density 
we also used the value of Kalberla et al. (2005), which are 
3.49, 3.41, and 3.43 x 10 20 cm" 2 for A1555, A1558 and 
SUZAKU J1231-1307, respectively. Since the statistics is 
limited we could not determine the metal abundance by 
spectral fitting, and therefore we fixed it at 0.3 Z Q which 
is the average value of the clusters in the distant cluster 
catalog compiled by Ota & Mitsuda (2004). We adopted 
the solar abundance table of Anders & Grevesse (1989). 
The redshift was fixed at z = 0.127 for A1555 and z = 
0.116 for A1558, respectively (Colafrancesco 2002). For 
the SUZAKU J1231-1307 we fixed the redshift at 0.1, since 
fitting is unstable when redshift is set to free due to low 
statistics. Then the free parameters are the normaliza- 
tions of the thermal and power-law components, temper- 
ature, and photon index. The results of spectral fittings 
are summarized in table 5, and the spectra are shown in 
figure 3, 4, 5, respectively. 

As a result, all the three models are consistent with the 
data of all the three sources in terms of the value of x 2 , 
except the power-law model against A1558 (5% probabil- 
ity of getting such a large x 2 )- However, in the case of 
A1555 the fit significantly improves when we fit by the 
thermal plus power-law model (x 2 /d.o.f = 27.8/44), com- 



pared with the pure thermal model (42.4/46). The F-test 
probability is 9.3 x 10~ 5 , indicating that the former is sig- 
nificantly favored than the latter. We examined the sensi- 
tivity of this result to the CXB fluctuation, by artificially 
increasing the CXB flux by 14.8%, which is la of the fluc- 
tuation (see section 2 for detail). The resultant F-test 
probability only slightly increases to 4.2 x 10~ 4 . We will 
discuss implications of those results in §4.1. 

In contrast, the spectrum of A1558 was not well fit- 
ted by the power law model, and x 2 value did not im- 
prove when we fit by the thermal plus power-law model, 
compared with the pure thermal model. Therefore A1558 
seems to be a well relaxed cluster. 

The spectrum of unknown diffuse source SUZAKU 
J1231-1307 is well fitted with all models. We will discuss 
this object later (see section 4.1). 

3.3. Total Mass estimation of A 1555 

Here we present a detailed estimation of the total mass 
of A1555, Mtotj because this source has been detected 
with the highest significance and the whole cluster is in- 
cluded in the field of view. 

In the isothermal /3-model(Cavaliere & Fusco-Femiano 
1976), the radial profile of X-ray surface brightness is writ- 
ten as 



S x (b) = S 



1 



-3/3+1/2 



(1) 



where 6, So, r c and /3 are the projected radius, central sur- 
face brightness, core-radius and the outer slope of density 
profile, respectively. 

We simulated the thermal emission from A1555 which 
follows (3 model profile using the ftool "xissim" , and com- 
pare it with the observation to estimate the value of /3 
and r c . As a spectral model, we used a single temper- 
ature thermal model with the parameters listed in table 
5. 

In figure 6, we show the observed surface brightness 
profile of A1555 with the results of simulations. First, we 
simulated with j3 = 0.5 and r c = 0.16 Mpc (shown as "single 
/3 model" in figure 6). This value of r c is an average of 
the clusters located at z = 0. 1-0. 15 in the distant cluster 
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SUZAKU J 123 1-1307 (PL) 



1 2 

Energy (keV) 



SUZAKU J1231-1307 (kT) 



1 2 

Energy (keV) 



SUZAKU J1231-1307 (kT+PL) 




1 2 

Energy (keV) 



Fig. 3. Spectrum of diffuse source, SUZAKU J1231-1307 with the power law model, thermal model, and thermal + power law 
model (left to right). Data points of XIS-FI are represented by black, while X1S-B1 by red. Solid line represents best-fit model. See 
Table 5 for fitted model parameters. 



ABELL 1555 (PL) 



ABELL 1555 (kT) 



ABELL 1555 (kT+PL) 



,0.01 



E10- 1 ? 



H 




4 



1 2 

Energy (keV) 



Energy (keV) 



1 2 5 

Energy (keV) 



Fig. 4. The same as figure 3, but for Abelll555. 
ABELL 1558 (PL) ABELL 1558 (kT) 



ABELL 1558 (kT+PL) 



0.01 



~10" 3 , 



10 




ppf— 



OS 



1 2 

Energy (keV) 



1 2 

Energy (keV) 



1 2 5 

Energy (keV) 



Fig. 5. The same as figure 3, but for Abell 1558. 
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Table 4. The best fit parameters of the CXB model for background spectrum. All errors are given at the 90% confidence limit. 



7 



^cooi f iVhot 1 fcThot 11 X 2 /d.o.f 



"PL 



Region 1 5^f£! 1.2!^ 4.7±t''( 0.26i^^ 110.1/96 

Region 2 5.8±g;| 1.3±g;| 0.29t°^ 93.6/121 

Region 3 5.2±g j 1.4±gJ 6.31^ 0.28^ : ^ 137.7/116 

Region 4 5.1±g'| 1.5±g'| bA±\\ 0.29±g g| 120.5/105 



* Surface brightness of power law component in 2.0 - 10.0 keV (in 
10 -8 erg cm -2 s _1 str -1 ). 

t Normalization in the thermal model, for the cooler and hotter galactic 
thermal components (in 10 — 2 cm -5 and 10~ 4 cm~ 5 , respectively). 
II Temperature of the hotter galactic thermal component (in keV). 



Table 5. The best fit parameters of the thermal, power-law (PL), and thermal plus PL model for the three diffuse sources: SUZAKU 
J1231-1307, A1555 and A1558. All errors are at the 90% confidence limit. 



Object 


Model 


Nth* 


kT T 




r s 


X 2 /d.o.f 


SUZAKU J1231-1307 


Thermal 
PL 

Thcrmal+PL 


0.9+83 
i-ol8:5 


3.7li.i 

°- J -5.6 
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catalog (Ota and Mitsuda 2004). The dispersion of r c is 
about factor two, and we confirmed that this uncertainty 
of r c does not significantly affect the estimation of the 
total cluster mass (about factor 0.5). The value of (3 is 
selected to fit the observation. This value is in agreement 
with the average of the distant clusters (Ota and Mitsuda 
2004) . It can be seen that the simulation well reproduces 
the observed surface brightness profile at the region except 
the central two data points. 

The central excess is often represented by an additional 
small scale /3-model component (double-/3 model; Jones & 
Forman 1984). We therefore performed a simulation by 
adding a small scale component of r c = 10 kpc//i and (3 
= 0.7, in addition to the single /3 model used above, as a 
test of the double (3 model. Those values of r c and (3 are 
averaged values of the small scale components in Mulchaey 
and Zabludoff (1998). The simulation result is shown in 
figure 6. The double f3 model gives a better fit, but is still 
discrepant with the innermost data by about two sigma. 
This may be just a statistical fluctuation, or may indicate 
a contribution from a AGN. Since it is difficult to further 
discuss the origin of this excess due to the low statistics 
of our data and limitation of spatial resolution of Suzaku, 
we conservatively consider a systematic error of a factor 
of about 2 for f3, from the dispersion of f3 in the distant 



cluster sample of Ota & Mitsuda (2004). 

Under the isothermal j3 model and also assuming that 
the cluster gas is in hydrostatic equilibrium, we can es- 
timate the total cluster mass from the value of r c , (3, 
and gas temperature (Cavaliere & Fusco-Femiano 1976). 
Using the single (3 model with f3 = 0.5 and r c = 0.16 Mpc, 
Mtot(< r 500 ) is estimated to be 7.43t2:2l5 4 o 5 x 10 13 M Q 
and r 50 o = O^lg ijgig^ Mpc, where the r 50 o is a radius 
where the ratio of total cluster mass density to critical 
density became 500. In the above mass estimation we used 
the value of gas temperature as 2.0i ) '33 keV which is de- 
rived from the single thermal component fit (see table 5). 
If we used the temperature derived from the thermal+PL 
fit, 1. 1 jlo'll keV, and the same value of (3 and r c assuming 
that thermal and non-thermal component have the same 
radial profile, M tot (< r 500 ) becomes 2.8! I'l-tl x 10 13 M o . 
The former errors in M to t and rsoo are statistical, simply 
estimated from statistical errors of kT and j3, while the 
latter ones are systematic, estimated from the systematic 
error of (3. The value of total cluster mass and tempera- 
ture are consistent with the known correlation between the 
total cluster mass and temperature (e.g., Ota and Mitsuda 
2004). 

If we adopt the double j3 model, the total cluster 
mass becomes 7.43 x 10 13 M Q (with thermal model) or 
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Fig. 6. The surface brightness profile of Abell 1555. The 
crosses are XIS-BI data. The dashed line represents best-fit 
single model. The solid line denotes double model while 
the dotted line is small scale component double-/? model. The 
typical background level has shown by dot-dashed line. 

2.8 x 10 13 M Q (with thermal + PL model). The differ- 
ence between the mass estimated with single j3 model and 
double (3 model is less than 1%. 

In the above mass estimation, we assumed that A1555 
is isothermal. We checked the effect of this assumption on 
the mass estimation by comparing the total cluster mass 
of several clusters estimated by Fukazawa et al. (2004) 
with the total cluster mass estimated by Vikhlinin et al. 
(2006), where the former assumed isothermal while the 
latter did not. These two differ by ~ 20%, indicating a 
similar level of systematic errors for our mass estimates. 

We also derive a luminosity of A1555 in 0.2 - 10.0 keV 
band from spectral fitting with thermal model. It is esti- 
mated to be 8.8 ± 1.6 x 10 43 erg s -1 . 

4. Discussion 

4-1- On the Evidence for Non-thermal Emission from 
A 1555 

Here we discuss the origin of the non-thermal emission 
implied for A1555. First we discuss the possibility that it 
is due to the IC scattering of CMB photons by accelerated 
electrons. From the T = 2.2 power-law component of the 
thermal+PL fit for A1555, we estimate the total energy 
of non-thermal electrons to be 3.1^'^ x 10 59 er g assuming 
z = 0.127, in the electron energy range 0.4-1.5 GeV cor- 
responding to the X-ray energy range of 0.5-7.0 keV. The 
total thermal energy of baryon gas when the structure for- 
mation occur is given by E haxyon ~ (3/4)(n B /^lo)M tot V^, 
where V c is circular velocity of the halo, and which can 
be calculated from spherical collapse model (Peebles 1980: 
Kitayama and Suto 1996) which is widely used in study 
of structure formation. If we use the total mass of A1555 
estimated with kT = 1.1 keV from the thermal+PL fit 
(see section 3.3) £ b aryon becomes 7 .OtJ^t^g x 10 60 erg, 
where the former error is statistic while the latter is sys- 
tematic which are estimated from the uncertainty of total 
mass (see section 3.3 for detail). Therefore the fraction of 



non-thermal electron energy becomes £ e ~ 0.044, which is 
consistent with the typical value discussed in the litera- 
ture (e.g., Fusco-Femiano et al. 1999, but see also Rossetti 
& Molendi 2004). 

From the flux of non-thermal PL component and the 
strength of magnetic field, B, we can estimate the flux 
of possible synchrotron emission due to the same elec- 
tron population which emits the non-thermal IC emis- 
sion. As already noted in section 3.2.1, A1555 has no 
radio counter part in NVSS catalog. Flux limit of NVSS 
is 2.5 mjy/beam and beam size is 45". Assuming that 
the synchrotron emission has the same spectral index with 
the PL component, T = 2.2, and extended about 1 Mpc, 
the expected brightness at 1.4 GHz is estimated to be 
2.2 mjy/beam with B = 1 /j,G. The typical value of B is 
0.1-1 fiG, and therefore the scenario of non-thermal IC 
emission is not inconsistent with the radio observation. 

We also discuss the detectability of the non-thermal 
IC emission with Fermi Large Area Telescope (Atwood 
et al. 2009). Theoretically it is reasonable to expect 
that the IC emission extends to the GeV energy band 
with a photon index of T ~ 2 (e.g., Totani & Kitayama 
2000). Extrapolating from the observed X-ray flux, the 
expected flux in the energy higher than 100 MeV became 
-FmoomcV = 3.9 x 10" 10 and 4.4 x 10" 11 photons cm" 2 s" 1 
for r = 2.0 and 2.2, respectively. These fluxes are about 
one order of magnitude lower than the sensitivity of Fermi, 
and the Fermi observation does not give a strong con- 
straint. 

It should also be noted that the flux of power law com- 
ponent in 13.0-40.0 keV is consistent with the flux upper 
limit obtained by HXD analysis. 

It has been known that relaxed clusters do have a cool 
core and a radially decreasing temperature at large radii, 
and merging clusters have shock heated structures. To 
approximate these effects on our results we fitted the two 
temperature thermal model to the spectrum of A1555. 
The fitting results are summarized in table 6. The fit also 
significantly improves comparing with the single temper- 
ature fit. The thermal+PL model has a smaller \ 2 value 
than the two temperature model, but the likelihood ratio 
is just a factor of about four. Therefore the two tempera- 
ture model provides another good description of the data, 
as an alternative to the thermal+PL model. The two tem- 
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peratures are determined to be 1.0_q 28 keV and 3.6 
keV, which are similar to the result of two temperature 
analysis of the Virgo and Centaurus cluster (Fukazawa et 
al. 2000). 

We also check a possibility that the origin of the power- 
law component is the background (or foreground) AGNs. 
We found two radio point sources at the outer region of 
A1555, NVSS J1232-1323 and NVSS 1232-1321 (Condon 
et al. 1998). Their flux density at 1.4 GHz are 7.0 + 0.5 
(NVSS J1232-1323) and 4.0+0.6 (NVSS J1232-1321) mJy. 
If we extrapolate their flux to the X-ray regime using the 
spectral energy distribution model of AGNs constructed 
by Elvis et al. (1994), the expected flux become 0.8 x 
10~ 14 - 1.0 x 10~ 12 erg s" 1 cur 2 in 2.0-10.0 keV. The 
estimated flux of PL component in the spectrum of A1555 



No. ] 

is 1.0 x 10~ 13 ergs -1 cnr 2 in 2.0-10.0 keV, and therefore 
the X-ray flux extrapolated from the flux of radio sources 
are consistent with the flux of PL component within the 
uncertainty. In conclusion wc can not reject the scenario 
that a AGN explains the non-thermal component in the 
spectrum of A1555. 

4.2. On SUZAKU J1231-1307 

We also detected interesting diffuse object, SUZAKU 
J1231-1307, which locates along the large scale filamen- 
tary structure. The spectrum of SUZAKU J1231-1307 
is well fitted by power-law, thermal, and thermal + PL 
model with similar X 2 value. If SUZAKU J1231-1307 is a 
galaxy cluster, the temperature is 3.7 keV and the total 
cluster mass is estimated to be 2.6l 2 Q x 10 14 M Q . Here 
we used the correlation between the total mass and tem- 
perature obtained by Vikhlinin et al. (2009). If the PL 
component really exists, the total energy of non-thermal 
electron becomes 5.0l 23 6 x 10 58 erg, which is about 15% 
of A1555. Further observation is required to determine 
the origin of this mysterious object. 

4-3. IC emission upper-limit associated with the large- 
scale structure as a whole 

If the non-thermal IC emission is emitted from the 
whole region of the large scale filamentary structure and 
extended more than XIS FoV, the flux of this emission 
has some contribution to the background flux. From the 
error estimation of CXB flux, we derive an upper limit to 
this component. As already noted in section 2, the sys- 
tematic error of CXB flux in one XIS FoV is estimated 
to be 9.65%. Combining this uncertainty with the sta- 
tistical error in CXB estimation (see table 4), total un- 
certainty of the CXB emission from all XIS FoV becomes 
5.4%. From this uncertainty and the average value of CXB 
surface brightness, the upper limit on the surface bright- 
ness of non-thermal IC emission is estimated to be 2.9 
xl0~ 9 erg cnr 2 s" 1 str" 1 . The total FoV of four XISs is 
1.1 x 10~ 4 str, and therefore the flux upper limit of non- 
thermal IC emission becomes 3.2xl0~ 13 erg cm" 2 s" 1 . 
This limit is in the same order of magnitude as the 
flux of power-law component in the spectrum of A1555, 
1.0xl0 _13 erg cm" 2 s _1 . 

5. Summary 

We performed the deepest X-ray observation of the su- 
pcrclustcr region including A1555 and A1558 with Suzaku. 
This region is interesting since large scale filamentary 
structures connected with the two galaxy clusters exist, 
but only RASS data was available in the X-ray band. 

We have detected the two galaxy clusters A1555 and 
A1558 for the first time in the X-ray band. Six point 
sources and one diffuse source (SUZAKU J1231-1307) 
were also detected. The total mass of A1555 was de- 
rived using the beta model. We examined the non- 
thermal nature of SUZAKU J1231-1307, A1555, and 
A1558, and there are no evidence for non-thermal emission 
for SUZAKU J1231-1307 and A1558. However, we found 
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that the A1555 spectrum cannot be fit by one-temperature 
thermal model. A thermal plus power-law model signifi- 
cantly improves the fit, and if this power-law component 
is the IC emission of CMB photons, the non-thermal elec- 
tron energy is estimated to be 2% of the thermal energy of 
the cluster, which is similar to the values often discussed 
in the literature for galaxy clusters. 

However, the A1555 spectrum can also be fit by a two- 
temperature thermal model, and the possibility of under- 
lying AGNs cannot be excluded. More detailed and sen- 
sitive observational studies are desired for better under- 
standing of this interesting region. 
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Table 6. Best fit parameters of two-temperature thermal model fitted to the A1555. All errors are at the 90% confidence limit. 
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